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Abstract 

The magnetostatic and magnetocapillary instability problems of isothermal incompressible and 
inviscid non-conducting liquid jets in a uniform magnetic field, is considered. The equivalence 
between static and dynamic approaches at the onset of instability and cut-off wavelength is 
shown. It is established that in the absence of electric currents the stability of permeable jets 
can be changed by the magnetic field. A new dispersion relation for magnetocapillary instability 
in such jets is derived. This relation differs from that given by Chandrasekhar. The existence 
of critical magnetic field which stabilizes jets with finite susceptibility is established. It is shown 
that the jet is stabilized by the field irrespective of its being para- or diamagnetic, but the extent 
of stabilization is different. 
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1 Introduction 



Joseph Antoine Ferdinand Plateau (1801-1883) was a Belgian physicist who is best remembered in 
mathematics for Plateau problem. He wrote his seminal book [1] being already blind for the last 
40 years of his life. The Plateau problem, in its brief formulation, is to find a surface of minimum 
area S given its boundary dQ. The problem relates to the principle of minimum free energy at 
equilibrium. In this context it may be required to find a surface with isotropic tension a which 
provides a minimum surface free energy aS. Rayleigh [2] gave a theoretical explanation for the 
instability of liquid cylinders that are longer than their circumference. Further generalization is 
called for if the excess free energy W of the cylinder comprises different types of energy that reflect a 
more complex structure of the liquid (e.g. elasticity [3] etc) and well as its capacity to interact with 
external fields. The threshold of static instability, which is associated with the change of sign of W, 
when the disturbed configuration becomes more preferable, is defined by 

W(kR,p) = -► k s = k s (p), (1) 

where A; is a wave number assigned to a small surface disturbance, R is a radius of cylinder, and the 
parameter p stands for the effect of an external interaction, which contributes to the jet evolution, 
e.g. fluid polarization in the presence of external fields. Note that the static cut-off wave number k s 
is related to the static cut-off wavelength A s by A s (p) = 2n/k s . 

Plateau instability being a static problem has also a dynamic aspect. Consider Rayleigh insta- 
bility in liquid jets of radius R and its corresponding dispersion equation s = s(kR,p). The latter 
determines the exponential evolution in time 0(r, t) ~ e st of all hydrodynamic functions 0(r, t) of 
the jet, with the growth rate s. Rayleigh's theory of capillary instability in liquid jets states that the 
maximum of the dispersion function s(kR,p) which corresponds to the wave number k max {p), gives 
rise to evolution of the largest capillary instability. The range of wave numbers k which contribute 
to the evolution of the instability is given by s > 0. Thus, the threshold of instability follows as 

s(kR,p)=0 - k d = k d (p). (2) 

Note that the dynamic cut-off wave number k d is related to the dynamic cut-off wavelength A d by 
A d (p) = 2-njk d . The value A d (0), which corresponds to free jet evolution is equal 2irR according 
to Rayleigh [2] . If there exist a critical parameter p = p cr such that for all p > p cr the expression 
s(kR,p) < holds, then for p > p cr the jet is stable for all wavelengths. The latter means that the 
liquid jet preserves its initial shape which is unaffected by small perturbations irrespective to the jet 
velocity. Thus, this conclusion must hold also in the limiting case of motionless fluid, i.e. 

A d (p)^A s (p). (3) 
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This identity reflects a deep equivalence between the static approach (the threshold of static insta- 
bility concerned with an excess free energy) and the dynamic approach (the bifurcation of the first 
non-trivial steady state of the inviscid hydrodynamic system). Thus, we come to Plateau problem 
for liquid cylinder with free surface in an external field. 

The capillary instability of magnetizable liquid jets in the presence of an axial magnetic field 
constitutes a classical example. Chandrasekhar [4] treated the hydrodynamics of magnetic liquid 
jets and derived their dispersion relation. This relation reads particularly simple for superconducting 
s sc {w) and non-conducting s nc {w) liquids (see [4], p. 545, formula (165) and p. 549, formula (205)) 



a, p stand for magnetic susceptibility, isotropic surface tension and density of the liquid, respectively, 
and po denotes the permeability of free space. I m {x) and K m (x) are the modified Bessel functions of 
order m of the 1st and 2nd kind respectively. In the superconducting limit the relation (4) discloses 
the existence of critical magnetic field H cr = H s /y/2 beyond which the jet is stable. In the non- 
conducting limit, this dispersion relation concides with that of Rayleigh [2] and 'in this limit the 
magnetic field has no effect on the capillary instability as should, indeed, be the case' (quotation from 
[4], Chap. 12, §112, p. 549). 

The latter result seems rather unexpected. Indeed, the absence of the electric currents in the 
presence of static magnetic field does not exclude coupling between the magnetic dipole moments of 
the liquid and the external field. Since the linear problem sets no constraints on the initial velocity of 
the jet, the above observation must be also correct in statics (Plateau problem) when the magnetic 
cylinder is destabilized. We proceed to verify our observation, first for the case of Plateau instability 
(see Sections 2, 3, 4) and then move on to the dynamic case. In Section 5 we give an accurate 
solution of the magneto capillary instability problem of isothermal, incompressible, inviscid and non- 
conducting jets, in the presence of a uniform magnetic field. We show that the dispersion relation 
differs from the one (4) found by Chandrasekhar [4]. The new relation accounts for the effect of 
magnetic fields and conforms with the solution of Plateau instability for magnetizable liquids. The 
reason for the discreapency in Chandrasekhar treatment arises due to incorrect boundary conditions 

which he applied to the magnetic field 1 . 

1 We quote from [4] ( Chap. 12, §112, p. 547) : 1 the field h must be continuous across the boundary '. The correct 
boundary conditions require continuity of the tangential component of the magnetic field and of the normal component 
of the magnetic induction (see equation (17) in Section 3 and equation (48) in Section 5.1 of the present paper). 




where w = kR and H s = a / '/io(l + x)R i s designated by Chandrasekhar as characteristic field, x 
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2 Free energy of liquid cylinder in the presence of magnetic field. 



The Plateau problem of static instability of a non-conducting liquid cylinder which is subjected to 
a uniform magnetic field appears, at first glance, to be a simple generalization of its counterpart 
in the absence of external fields. However, deeper consideration shows that the presence of the 
field complicates considerably the physical picture and computational procedure. A fundamenatal 
question arises concerning the correct definition of the excess free energy W which must be minimized 
via variation of the cylinder shape. 

Consider an isothermal liquid cylinder in a uniform magnetic field Ho that is applied in free 
space along the cylinder axis. The magnetic susceptibility \ °f the cylinder is assumed isotropic, 
independent of magnetic field, and satisfying the thermodynamic condition % > — 1 [5]. When the 
liquid cylinder is undisturbed the total free energy F° of the system is given by 



F» = £»-±^. dv, (5) 

cyl 



where the integral represents the volume ttR 2 L, enclosed by the area dQ,® yl which is occupied by 
undisturbed cylinder. The term £® = a j dn0 ds = 2-noRL stands for the surface free energy of the 
undisturbed cylinder, where R, L and a denote its radius, length and surface tension respectively. 
Deformation of the cylinder shape changes the field H(r) over all space M 3 , i.e. in both the internal 
domain Q cy i and its complement (the exterior domain) M 3 \ £l cy i. Following Plateau, we assume 
conservation of the cylinder volume 

f dv= f dv . (6) 

^^cyl J^cyl 

The total free energy F of the disturbed cylinder takes the following form, 



F = S s - (1 + X )ff ■ I -'"U\r)dv -ff.f ~H\r)dv + f ■ / H*dv 



(7) 



\yl Z ^ 3 \Ji CB i 

where m H(r) and ex H(r) denote the internal and external magnetic fields, and £ s = a f dn ; ds. The 
excess free energy W of the system is defined as, 

W = F — F° . (8) 

From the mathematical standpoint, the variational problem for minimization of W, supplemented 
with constraint (6) for all smooth surfaces dQ cy i, is known as the isoperimetric problem. 

The cylinder instability can be studied assuming small perturbation in its shape. In this case the 
Plateau problem becomes solvable in closed form. 
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Let the extent of deformation be characterized by a length (o, such that (o/R = e 1. Then 
the following approximation holds 

/ ds - [ ds ~ Ae 2 RL , (9) 

Jan cyl Jan« yt 

where A is a dimensionless parameter dependent on the deformed geometry. 

Here, the fields m H(r) and ex H(r) which must satisfy Maxwell equations can be represented as 
small perturbations of Ho, 

in H(r) = H + in H 1 (r) = (H + m H l z , ' m H l r ) , ex H(r) = H + ex H 1 (r) = (H + ex #i, ex H l r ) , (10) 

where according to the assumption eCl the following approximations apply (see Setion 3) 

{ m Hl ]n Hl ~Hl) = {'"hi, in hl, "hi} x e X H . (11) 

The dimensionless fields m,ex hl z (r, z) are dependent on the coordinates. Inserting (11) into (7), and 
performing the integration, we evaluate W given in (8) as, 

W = Aae 2 RL - ^-U , (12) 

where use was made of (5) and U is given by (see Appendix A) 

U = (1 + X) / { ( n Hl) 2 + <!'>H}) 2 }dv+ ! { {~Hl) 2 + {~H}) 2 }dv + 
Jn cyl L J JM.»\n cyl L J 

2H ((l + x)[ m Hldv+ I ) . (13) 



Q C yl ^R 3 \n c yi 

We specify the commonly used harmonic deformation of the cylinder as r(z) = R + Co cos kz, where 
k = 27r/A, A being the disturbance wavelength. By virtue of translational invariance of the problem 

in ' ex /i^(r, z + A) = > n ' ex hl. z (r, z) (14) 

we set L = A and evaluate the free energy per unit wave length. For this type of deformation the 
parameter A is known [1], 
r 

a ds- 2naRL = a-^§L (w 2 - l) — ► A = - (w 2 - l) , w = kR . (15) 

Jdn cyl 2ti 2 

Next we proceed to solve the boundary problem and get the distribution of the magnetic fields. 



3 Boundary problem and its solution 

The magnetostatics of the disturbed liquid cylinder is governed by Maxwell equations for the internal 
in H(r) and external ex H(r) magnetic fields 

rot in H = rot ex H = , div in B = div ex B = , (16) 



where m B = + x) m H and ex B = /io ex H denote internal and external magnetic inductions, 

respectively. Equations (16) must be supplemented with boundary conditions (BC) at the interface 
r = R, 

< in H, t) = < ex H, t) , < in B, e) = < ex B, e) , (e, t) = , (17) 

where t and e stand for tangential and normal unit vectors to the surface, respectively. Since the 
surface deformation is small, linearization can be applied, 

t r = -e z = d(/dz , t z = e r = ^l-(d(/dz) 2 ~1. (18) 

A standard way to solve the problem is to introduce the magnetic potentials <£in(r) and ( I > ex(r) which 

are defined as '"H^r) = -V$i n , ex H 1 (r) = -V$ e x, where |V3>j n |, |V$ex| < H . These potentials 

satisfy the first two equations in (16). The last two equations in (16) yield, 

B 2 (b- 8 2 (b B 2 1 8 

^ + A 2 $ in = 0, ^ + A 2 d> ex = 0, A 2 = ^ + if , (19) 

oz az ar^ r or 

where A 2 is the two-dimensional Laplacian. Reformulation of the BC (17) for & ex (r), &- m (r) gives, 

d$ ex d$ m 8<s> ex d<s> m d( 

-^z- = ^T> -dr--( 1+X) -dr-= XH d-z- (20) 
Using $>\ n (r, z) = 4>\ n {r) sin kz and $ex(^ 5 z) = 4> ex (r) sin fez we find 

(A 2 - k 2 ) ex = , (A 2 - k 2 ) in = (21) 

with BC at r = R 

*« = *n, (l + x)^~^ = xHkCo. (22) 
The solutions of (21), which satisfy BC (22) and are finite at r = and r = oo, are obtained as, 

4>Ur) = CoXH b -^^Io(kr), 0«(r) = CoxH ^^-K (kr) , (23) 
w w 



where 



2 K (w) 2 h{w) 



K^i x) = ™ 7f ( r , c ( ro > X) = ™ 7^7 r ' T ( ro ' X) = 1 + X^h(^)K (w) . (24) 

Recalling the definition (11) of the dimensionless fields m,ex h}. z (r, z), the following final expressions 
are obtained, 

'"hi = — b{w, x)Io(kr) cos kz , '"hi = — b(zu, x)h(kr) sin kz , 

ex hl = -c(zu,x)K (kr) cos kz , ex hl = c(w, x)Ki(kr) sin kz . (25) 

Bearing in mind that the function Q(w) = tuIi(tu)Kq(w) is monotone growing, at the positive half 
axis and is bounded < Q(w) < 1/2, we conclude that the fields m,ex /i^ z (r, z) are free of singularities 
in the thermodynamically relevant region x > — 1. 



4 Plateau instability 

In this Section we calculate via equation (11) the excess free energy W in terms of the dimensionless 
fields ex hl r , m h\ r . Inserting (25) into (13) we get (see Appendix A) 

W = — - — e ■f(w,x,H ), f(w,x,H ) =w -1 + XOT 7f7 r~ , (26) 

where the dimensionless excess free energy f(zu, x, Hq) can be defined by introducing the character- 
istic field Hch 



In(w)Kn(m) 



H, 



Ch 



a 



(27) 



Notice that Hch differs from H s introduced by Chandrasekhar in (4). Formula (27) serves for both 
paramagnetic (x > 0) and diamagnetic (x < 0) liquids, since the latter case does not lead to an 
imaginary expression in (27) due to the term x 2 in (26). 

Consider the field H cr (w,x) satisfying f(w,x>H cr ) = and call it critical field. Beyond this 
field, Hq > H cr , the cylinder is stable. The expession for the critical field is 



Hcr(™,X) 



H C h l-w 2 T(w, x) 



zu 



(28) 



X Iq{w)Kq{w) ' 

Figures 1, 2 show plots of h cr = H cr /Hch for strong (x>l) and weak (|x| <C 1) magnetic liquids. 

strong paramagnetism 
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Figure 1: A plot of h cr versus k s R for large positive susceptibilities : x = 1> 10) 10 2 , 10 3 , 10 6 , from 
right to left, respectively. 

The corresponding asymptotics for H cr (vu,x), in the case of weak x *C 1 and strong x 3> 1 
magnetic susceptibilities and in a long wave limit w —* 0, which are derived in Appendix B, are 
presented below 
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Figure 2: A plot of h cr versus k s R for small positive and negative susceptibilities: x = 
±10 -4 , ±10 -3 , ±10 -2 , ±10 , from right to left, respectively. The plots corresponding to diamagnetic 
and paramagnetic cases nearly coincide. 



To the end of this Section we present physical arguments which justify our conclusion regarding 
stabilization of the permeable jets for both paramagnetic and diamagnetic liguids. The orientation 
of the dipole moments along the applied magnetic field in the paramagnetic liquid enforces the 
rigidity of the jet under small disturbances. In the case of diamagnetic liquid the increase of the 
jets' stabilization is due to orientation of the dipole moments in the opposite direction which also 
enforces the rigidity of the jet. The similar increase of the rigidity, and its corresponding influence 
on the stabilization, was found recently [3] in the liquid crystalline jet due to elasticity of the media. 
Although the critical fields for paramagnetic H cr (w,x) and diamagnetic H cr (w, — x) liquids are 
different the following universal relation holds, 



5 Hydrodynamics of Non— Conducting Jet in a Magnetic Field 

Consider an isothermal, incompressible, inviscid and non-conducting jet in the presence of a magnetic 
field Ho applied along its z-axis. The deviation from initial values of the pressure is defined as 
P{ n = P m — P n , where P n is the unperturbed pressure within the cylindrical jet. The deviations 
of the internal and external magnetic fields are defined as '"H 1 — n H — Ho and ^H 1 = ex H — Ho, 
respectively. The governing equations of magnetohydrodynamics which are given by 



where B 1 = l/y/I (l)K (l) ~ 1.3697, B 2 = h{l)K (l) ~ 0.2379. 



H 2 cr (w, X ) Hi 





(32) 



divV" = 



div a U = 



P 




(33) 



(34) 
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can be simplified considerably by applying Maxwell equations (34) for non-conducting media (71,, = 
7ex = 0) to Navier-Stokes equation (33). y m and 7 ex denote conductivities of the jet's interior and 
exterior, respectively, stands for magnetic stress tensor, and V m is local fluid velocity. Finally 
the magnetohydrodynamic problem is decoupled into the hydrodynamic and magnetostatic parts : 

divV in = , ^— = -igrad P m , and div a H = , rot Q H = , a = in, ex . (35) 
at p 

A standard way to solve the boundary problem (35) is to introduce the Stokes stream function 
*(r, z, t) and magnetic potentials &- m (r, z, t), <£ ex (r, z, t) 

V r = ——, V z = -—, Q H = Ho + a H , "H^-V^, a = in, ex. (36) 
r oz r or 

This gives the following governing equations 

1 dP m 1 a 2 * 1 dP m 1 a 2 * / A d 2 \ a 

+ 71^ = 0. T^--r^ = 0' A 2 + ^ *« = , a = .n,ex. (37) 



p dz r drdt ' p dr r dzdt ' \ dz 2 1 

These equations must be supplemented by four boundary conditions which are derived in the next 
Section. 

5.1 Boundary conditions 

It is neseccary to apply boundary conditions (17), which are imposed on H and B, as well as those for 
the hydrodynamic variables. First, the velocity V r must be compatible, at r = R, with the assumed 
form of the deformed boundary d(/dt. Second, at the free surface of a liquid jet the jump in stress 
must be balanced by Laplace pressure [5], 

[T rz ]ex e z + [ T rr\ex e r = 2aUe r , [T zz ]'" x e z + [T zr ]'" x e r = 2aHe z , (38) 
where [7)fc]^ = Tj£ — Tj^ and H is mean surface curvature decomposed as 

W = «b + «i, «o = ^. « 1 = -i(^ + g) Ke =|. (39) 
By virtue of (18) we get 

[T rr C - 2aH = [T rz t ^ , [T zr t = ~ 2a7i) g , or 

[T rr t - 2aH = ([T zz t - 2aH) (g) ' , [T-J^ = ([T„fi - 2aU) g , (40) 

where 

K n r = P m + ^{{ ' m H z f - ( ' m H r ) 2 } , TZ =P m + ^ {( ' m H r ) 2 - ( ]n H z ) 2 } , = -p >"H r >"H Z , 

T ex = ^ {( ex^ )2 _ ( e Xjffr) 2| ? j« = { ( ex^ _ ( ex^ )2 | _ j« = ex^ ex^ 



Recalling that d(/dz oc e and 



m H z -H = m Hl oc e , ex H z - H = ex ^ oc e , m H r = m H} oc e , ex H r = ex Hj oc e , 
( in tf z ) 2 = H 2 + 2tf in #i + 0(e 2 ) , ( ex H z ) 2 = H 2 + 2H ex H l z + 0(e 2 ) , (>"H r ) 2 = ( ex tf r ) 2 = 0{e 2 ) 



we get within the 1st order approximation in e 
[Trr]Z ~ 2(J ' H 



pi" _ 2 <jH + {|( in tf 2 ) 2 - ^( ex tf ^ 2 



2 _ M 

2 v J ' I 2 ' 2 

MOX T_r2 i tj /'..in ,. exT_j-l\ i /o^2\ 



H 



M0 | 

2 ' 



fl,)^ - '-i in //,.) 2 } 



pi" _ 2aW + ^# 2 + iJ (/* in ^i - Mo ex ^l) + 
[«-2oW = P i "-2^+{|( in F r ) 2 -^( ex J ff r ) 2 } + {^( ex J ff 2 ) 2 -|( in F z ) 2 } 

pin _ 2aH _ _ Hq ^ \n H l _ M + ^ 

[T zr C = ^ ex H r ^H z -^ m H r ' m H z = H Q {^^H l r -^ n Hl)+0(e 2 ). 



(41) 

(42) 
(43) 



Assuming P{" oc e and combining the first equation in (40) with (41) we obtain 

MX rr2 



pi" _ 2aH + q±H 2 + H (/x - // ex ^) = 



P Q " - 2aH + l^oxHi = , 
P{ n - 2<rWi + # (/* '"Hi - fi ex Hl) = 

This gives the unperturbed pressure within the cylindrical jet as, 



pin _ a ,, Mo#o 

Combining the second boundary condition in (40) with (42), gives 

{fio ex H r - y, ™H r ) H = (p™ - 2aH - ^H 2 ) 
Both equations (44) and (45) lead to the conclusion 



dz 



/j,q ex H r — jx m H r 



dz 



that coincides with the second static boundary conditions in (17). Thus we arrive at 

d( 



^ n + (J (^ + 0) +Fo ^ in ^ 1 - /i ° eX ^ =O 



V r = 



(1 + X ) m H r - ex H r = X H 



in ttI ex ttI 

g~ z , H z - H z , 



or using the notations of (36) 



dr 



(1 + X) 



d^ = xH d( 

dr dz 



d$ m 

dz dz 



ldV dC 

r dz dt 



(44) 



(45) 



(46) 



(47) 
(48) 



(49) 
(50) 
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Assuming that an axisymmetrical disturbance, characterized by a wavelength 2ir/k, increases expo- 
nentially in time with the growth rate s, gives, 

{$ in , $ ex , vl/, C, P[ n } = {i Mr), i 0ex(r), i Y>(r), c(r), p(r)} x e st+lkz . (51) 

Consequently, the following boundary conditions hold 

<9</>ex n . v ^in , „ , , 

-q^T - (! + X) = fc X#0<? , 0ex = 0in • 

The presence of the field-dependent term, in the first equation of (52), indicates that the coupling of 
hydrodynamics and magnetostatics in the boundary conditions must break the Rayleigh dispersion 
relation (4) for non-conducting jet. 

5.2 Dispersion relation 

Inserting (51) into (37) results in the following amplitude equations, 

^ + S pt=o, k P + sp-^ = 0, (A 2c - k 2 )^ = , (A 2c - k 2 )<p ex = . (53) 
k Or r r or 

The latter have fundamental solutions that are finite at r = and r = 00, as follows 

Y>(r) = Aikrhikr) , p(r) = -A lS pkI (kr) , </> in (r) = A 2 / (A;r) , (f) ex (r) = A 3 K (fcr) , (54) 

where Aj are three indeterminate coefficients. Substituting (54) into (52) we get the following 
dispersion relation 

2 a whjw) 2 2 3 p H 2 h(w)K {zu) 

S = D3 T ( \ (- 1 ~ ^ ) ~ X ^7 — ^ — r . (55) 

I (w) R 2 p T(zu,x) 

(see Figure 3). The last expression clearly differs from Chandrasekhar's result (4). Moreover, (55) 
can be recasted as 

R 6 p I (w) 

where f(w, x, Ho) can be recognized as the dimensionless excess free energy found in (27). Expression 
(56) confirms the equivalence (3) between the static and dynamic approaches in the problem of cut- 
off wavelength since both non-trivial zeroes, Wd and w s , of s nc (wd) and f(w s , %, Ho), respectively, 
coincide. 

It is noteworthy that in the limit \ — > 00, 

which coincides with the dispersion relation (4) for a superconducting jet in the longwave limit 
w -> 0. 
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Figure 3: A plot of the dispersion relation s(w s ) for positive susceptibilities x = 0, 1, 3, 10, from right 
to left. 

6 Conclusion 

We have revised the theory of capillary instability of isothermal incompressible inviscid and non- 
conducting liquid jets in a uniform magnetic field which was treated by Chandrasekhar [4] about 
50 years ago. The main result of [4] states that for non-conducting jets the field has no effect on 
their stability and the dispersion relation is given by Rayleigh's theory for free jets. This statement 
contradicts the free energy approach (Plateau problem) in the framework of static instability of a 
magnetazible liquid cylinder in a uniform field. 

We have found a new dispersion relation for magnetocapillary instability in such jets. The relation 
differs from that given by Chandrasekhar and agrees with the result for the cut-off wavelength 
obtained from the static approach. This reflects a deep equivalence between the static approach 
(the threshold of static instability which is concerned with an excess free energy) and the dynamic 
approach (the bifurcation of the first non-trivial steady state of the inviscid hydrodynamic system) . 
In this context, the existence of critical magnetic field which stabilizes jets with finite susceptibility 
is established. It is shown that the jet is stabilized by the field irrespective of its being para- or 
diamagnetic, but the extent of stabilization is different. 
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A Contribution of the magnetic field to free energy 



Evaluate the contribution hqU /2 of the magnetic field inside Q cy i and outside M 3 \f2 C2/ ; of the disturbed 
liquid cylinder to the excess free energy W (see equation (12)) 



U = X 



inTj2 



H z (r)dv 



cyl 



cyl 



cyl 



Hldv 



+ f ' ,n H 2 {r)dv+ f ex H 2 {r)dv- [ H$ 

Jn cyl J^\n„ d Jrz 

{ (H + in Hl) 2 + (<"Hl) 2 }dv- [ Hi 

cyl 

f {{HoV n Hl) 2 + ^Hl) 2 }dv+ f {{H Q+ ~Hl) 2 + ^Hl) 2 }dv- f H% 

(! dv- f dv)+ f {2H m H\ + ('"Hi) 2 + ( ,n H}) 2 } dv 
\Jn cyl Jn° yl J Jn cyl ^ 

( {2H m H\ + ("Hlf + ('"Hi) 2 } dv + / \2H ^H l z + ( ex ^) 2 + ( ex ^) 2 | dv 
Jn cyl ^ > JvL 3 \n cy i 1 } 

= (1 + X ) I |2tf '"Hi + ("Hlf + ( n Hl) 2 } dv+ f [2H, ex #i + ( ex ^) 2 + ( ex i^) 2 ) 
Jn cyl <- > JR 3 \n cy i 1 J 

= (1 + X ) f { (»Hl) 2 + ("Hl) 2 \dv+ [ { (*Hl) 2 + ^Hl) 2 \dv + 

ldv+ I 
Jr 



Hi 



dv 



'^cyl 

2H \(l + x) 
Hence, we conclude that 



ex ttI 



lM?\n cy i 



Hldv 



U = [C x e 2 X 2 + 2C 2 e X ) H 2 , 



(Al) 



where 



J_«~ = 1±A 
irL irL 



[ { ("hlf + (X) 2 } dv + -L / { ^h\f + (-hlf} dv , (A2) 

cyl \ cyl 



irL 



C 2 



cyl 

!+X / in, "J 



J 



h r dv + 



1 



exil 



jR3\ ncy[ 



(A3) 



Recall the expression (12) for W 

W = ^^ 2 -l)ae 2 RL-( Cl e 2 x 2 + 2C 2 ex)^ 



2 

nLaR 



J ^- 1 )-(^ + *S) 



C 2 \ XW)#o 



21 



<T R 



(A4) 



Making use of the integration 

/ [I$(v) + lf(v)] vdv = uh(u)h(u) , / [K 2 (v) + K 2 (v)] vdv = uK^K^u) , 

J Ju 

ru roo 

/ Io(v)vdv = ul\{u) , / Ko(v)vdv = uK\(u) , 

■/ Ju 
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and inserting (25) into (A2) and (A3) gives, 



-L/ o {("hl) 2 +(<X) 2 }dv = ^^.2n^f\l 2 (kr) + I 2 (kr)}(kr)d(k 

cyl " 



b^R 2 

w 2 

„2 j t>2 roo 

— 2vr^ / \K 2 (kr) + Kl{kr)} (kr)d(kr) 
ttL 2 zu z J kR 

i^-wKoWKW = ^R 2 



1 

7rL 



/ {r^) 2 +( e ^) 2 }d, = 

JR :i \n° cyl L J 



*,I (w)I 1 (w) = mR ^ TM J 



71 L Jn cyl 



dv = 



-I 



b R 2 f L rkr(z) 

• 2-7T — 7T / dz cos kz / Io(kr)(kr)d(kr) 

ttL w £ J J 

26 R 

/ cos kz ■ (R + Co c °s &z) • Ji (kR + &:Co cos fcz) dz 
Jo 



ex h\dv 



L w 



c R 
■ — • 2vr— 
ttL vj z 



2 pL f>oo 

dz cos kz / Ko(kr)(kr)d(kr) 

./fcr(z) 



2c -R f^ 1 

= — — — / cos kz ■ (R + Co cos kz) ■ K\ (kR + feCo cos kz) dz . 
L w J 

Expanding I\[x + e) and K\(x + e) in the vicinity of x = 0, up to the first order in e, gives 
I^x + £ ) = I^x) + (J (x) + J 2 (s)) I , + e) = K^x) - (K (x) + K 2 (x)) | . 



Finally we arrive at 



-/ 



26 -R /"^ 

-— — / dz cos fcz • <j RI\(w) + Co cos fez 



I\(w) + w 



I Q (w)+I 2 (w) 



\ f ex h]dv 
Jm?\Q, cyl 



26 R 

— Co 
L w 

2c R rL 
L w 

2c R 



h(vj) + w 



I (zu)+I 2 (w) 



I dz cos 2 kz- = —eb(w,x)R 2 Io(' u ^) 
Jo 



J dz cos kz ■ IrKi(w) + Co cos kz K\(w) — w 



K Q (w) + K 2 (w) 



L w 



Co 



K\ (vj) — w 



Kpjw) + K 2 (w) ] f L 2 2 

— — - / dz cos kz = ec(w,x)R Kq[w) . 

2 J jo 



Thus, we get 
1 



ttL 
1 



3 D 2 f / q(ct)Kq(ot) 

r(ro,x) 



(1 + X) 



+ 



I (tu) K (w) 



w K 



T(zu,x) 



—C 2 = -ew 2 R 2 

that consequently gives, 
irLaR 2 



(1 + X) 



w = 



I (zu)K (w) _ I (zu)K (zu) 
T(w, X ) T(w,x) 



Ci 



2d2 i o('^)Ko(zu) 



e 2 -f(w, X ), S{vj,x)=vj 2 -l- + 2 



r(u7,x) 

C 2 \ X/^o^o 



e7rL / cri? 



Insertion (A5) and (A6) in the latter expression, we obtain 



f(™,x) = w -1 + x ™ : 



2 2 ^oRHq I g {w)K q {w) 



v T(w,x) 



(A5) 
(A6) 

(A7) 

(A8) 



15 



B Asymptotics of expressions 

Consider the expressions obtained in Appendix A and evaluate its asymptotics. 

• w^O , \ X \ < 1. 

where 7 = 0.577216 is Euler' constant. The corresponding asymptotics for the dimensionless excess 
free energy f(w,x) and the critical field H cr (w, x) read 

f(w, x) = w 2 -l + x^ftw) (l - \ w 2 (3(w)) (J^ 2 , (Bl) 

H cr (zu, X ) = Hch * (l-| CT 2 ln OT ) , (B2) 

w y/— xlnro V 2 / 

• zu — ► , x ~^ 00 • 

/<„,X) ~^-l + 2 (l + ^)(^) 2 - = ■ (B3) 

• 07 — > 1 

/fox) - ^ 2 - 1 + : X^#7n f^) 2 ' H„(w t x) = B 1 H Ch jBr^VT=rf, (B4) 
1 + X-ti(lj-Ko(lj \ H ChJ V X 

where £1 = l/v%(l)^o(l) - 1-3697, 5 2 = Ji(l)K (l) ~ 0.2379. 
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